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Abstract 

We consider weak solutions for a diffuse interface model of two non-Newtonian 
viscous, incompressible fluids of power-law type in the case of different densities in 
a bounded, sufficiently smooth domain. This leads to a coupled system of a nonho- 
mogenouos generalized Navier-Stokes system and a Cahn-Hilliard equation. For the 
Cahn-Hilliard part a smooth free energy density and a constant, positive mobility is 
assumed. Using the L°°-truncation method we prove existence of weak solutions for a 
power-law exponent p > d = 2,3. 
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1 Introduction 

We consider a two-phase flow of two incompressible non-Newtonian fluids in a bounded 
domain. In the case of classical sharp interface models the interface, separating the two 
fluids, is modeled as a (sufficiently smooth) hypersurface. But these models do not allow to 
describe flows beyond the occurrence of topological singularities, e.g. when droplets collide 
or pinch off. In the following we will consider a diffuse interface model for such two-phase 
flows, where the fluids are assumed to be partly miscible and the sharp interface is replaced 
by a thin interfacial region, where a scalar order parameter p changes smoothly, but rapidely 
between two distinguished values, e.g. ±1, that describe the separate phases. More precisely 
we consider the following system, which couples a nonhomogeneous generalized Navier- 
Stokes system and a Cahn-Hilliard equation: 

gdt'v + {g-v + J)) • Vv — div S((^, Dv) -|- Vtt = —eo div (jVip (g) V:^), 

divv = 0, 


( 1 . 1 ) 

( 1 . 2 ) 


1 



(1.3) 

(1.4) 


dt(p + V • Vv9 = mAfi, 

H = e~^f{(p) - eoA(p 

in a space-time cylinder Qt = x (0,T), where J = —and H C d = 2,3, is a 
bounded domain with C^-boundary together with suitable boundary and initial conditions 
specified below. Here v: Qt —>• is the velocity of the mixture of the two non-Newtonian 
fluids, which is defined as the volume average of the individual fluid velocities, cf. [3], 
tt: Qt —>■ M is its pressure, Qt —)• M is an order parameter related to the volume 
fractions of the fluids (e.g. ip is the difference of the volume fractions of both fluids), 
/: M —>■ M is the homogeneous free energy density of the fluid mixture, p = p{ip) is the 
density of the mixture, depending explicitely on ip, and p: Qt —)• K is the chemical potential 
of the mixture. Moreover, S((/9, Dv) is the viscous part of the stress tensor due to friction 
in the fluid mixture, which will be specified below. Finally, Sq > 0 is a constant related to 
the thickness of the diffuse interface and m > 0 is a mobility coefficient, which is assumed 
to be constant as well. For simplicity we will assume eq = m = 1. But all results in the 
following remain valid for general eQ,m > 0. 

The model above is a non-Newtonian variant of the model derived in [4], where the 
constitutive assumption S((/j,Dv) = 2zx((^)Dv is made. A prototypical example for the 
following is 

(1.5) S((/9, Dv) = 2zx((^)|Dv|^“^Dv for some p > 1 

and a suitable positive M —>• M. Such power-law models for non-Newtonian fluids are very 
popular among rheologists The value for p can be specified by physical experiments. 

An extensive list for specific values for different fluids can be found in [9]. Apparently 
many interesting p-values lie in the interval [|,2]. The mathematical discussion of power- 
law models for non-Newtonian fluids started in the late sixties with the work of Lions and 
Ladyshenskaya (see [26l [271ES] and [28] ). A first systematic study can be found in [29| . 

In the case of (|1.5p the derivation of (ll.ip - (jl.4|) is precisely the same as in jH Section 2]. 
One just has to choose S((/?, Dv) as above, which guarantees the validity of the local dis¬ 
sipation inequality. Moreover, let us note that in the derivation of the model in [3] it is 
assumed that 

(1-6) = + 

where gi, Q 2 are the specific densities of the two (separate) fluids. Then §^(<p) is constant 
and g(ip) solves the continuity equation 

(1.7) 5t£»((/?)div(^>((^)vJ) = 0, 

where J is a flux of the fluid density due to diffusion relative to the flux pv caused by 
convection. Here (^ = (^2 ~ T’l is the difference of the volume fractions ^> 1,^2 of the 
fluids. Physically ip should only attain values in [—1,1], which guarantees that g(ip) is 
positive. Since the density is a function of the order parameter, the study share certain 
similar features with the analysis of quasi-compressible fluids, see for instance m and the 
references therein. We note that this diffuse interface model corresponds to a two-phase 
flow with a sharp interface T{t) separating two immiscible incompressible non-Newtonian 
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fluids. Here the surface T{t) gives rise to a surface energy given by a constant surface tension 
coefficient times the area of T(t). Note that no variable surface tension or curvature effects 
are taken into account in the surface energy. We refer to [H Section 4] for the relation 
to sharp interface models in the Newtonian case, which can be modified to the present 
situation. 

In the following we want to construct weak solutions of the system above for arbitrary large 
times 0 < T < oo. But, since a comparison principle for the fourth order Cahn-Hilliard 
system (ll.3l) - (ll.4l) is unknown and we will assume that /: M —)• M is a suitable smooth 
function, we are not able to prove that (p attains only values in [—1,1]. Let us note that 
in the case of Newtonian fluids (i.e., p = 2) the existence of weak solutions of the system 
above for large times was proven in laE]. In these contributions either /: [—1,1] —)• M is 
assumed to be singular at 9[—1,1] or the mobility m is a degenerate function of (p. In both 
cases one obtains that (p € [—1,1] almost everywhere and one can assume that (II.6p holds 
true. Since in our setting we are not able to show ip £ [—1,1] and p{ip) in (II.6p becomes 
negative outside of [—1,1] unless pi = P 2 , we cannot assume (II.6p for all values of ip and 
p{ip) has to be modified outside of [—1,1] such that it stays strictly positive. But then ^ 
is no longer constant for all values of ip and g solves only 


( 1 . 8 ) 


dtg + div(^)v + J) = i?, where 


R = -v|^ • V/r, 

dip 


instead of (|1.7I) . Here R is an additional source term, which vanishes in the interior of 
{ip G [—1,1]} if (II.6p holds for all ip G [—1,1]. In order to obtain a local dissipation 
inequality and global energy estimate the equation of linear momentum (11.11) has to be 
modified to 


(1.9) gdt'v + {g'v + J) • Vv + R— — div S{ip, Dv) + Vvr = —eo div (Vip (8) V(^). 
This modification guarantees that 


di 


£'|vl' 


V 


+div I (^3v + J) - ) —div(S((/?, Dv)v— ttv)- S((/?, Dv) : Dv = —Eq div(V99(8iV(/?)-v 


and the (global) energy identity 


dt 


[ ^ 

In 2 


■dx + J ^ _ J S{ip, Dv) : Dvdx — J mjV/r]^ dx 


for every sufficiently smooth solution of 

( 1 . 10 ) < 

where J = together with the boundary and initial conditions 


gcltv + (gv + J)) • Vv + — div S((/9, Dv) + Vtt = —eo div (Vy? (g) V(/?) , 

div V = 0, 

dtip + v-'Vip = mAp, 

= ^0 V'(¥^) - eoAip, 


( 1 . 11 ) 


vjan = 0 

djv^lan = dji/pldn = 0 
v(0,-) = Vo 
(^(0,-) = ipo 


on do, X (0, T), 
on dQ X (0, T), 
in Q, 
in Q. 
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Here AA denotes the exterior normal of dQ. Finally, we note that (jl.9l) is equivalent to 
(1.12) dtigv) + div(v 0 (gv + J)) — R— — div Dv) + Vvr = —eg div (Vip 0 V<y9) 

due to (11.81) . This reformulation will be used for the definition of weak solutions below. 

In the following we will prove existence of weak solutions for (jl.lOji - dl.lljl . So far 
(jl.iop - dl.lll) was only treated in the case that both fluids have the same densities, i.e., 
p{ip) = const.. First analytic results in this case were obtained by Kim, Consiglieri, and 
Rodrigues [3]. They proved existence of weak solutions if p > d = 2,3. In this 

case monotone operator techniques can be applied. In [2l] Grasselli and Prazak discussed 
the longtime behavior of solutions of the system in the case p > d = 2,3, assuming 

periodic boundary conditions and a regular free energy density. For the same p results on 
existence of weak solutions with a singular free energy density / and the longtime behavior 
were obtained by Bosia m in the case of a bounded domain in M^. Finally, existence of 
weak solutions was shown by Abels, Diening, and Terasawa [3] in the case that p > 
using the parabolic Lipschitz truncation method for divergence free vector helds developed 
by Breit, Diening, and Schwarzacher which is the same range for p as for a single 
power-law type fluid, cf. Diening, Ruzicka, and Wolf [17] (the same bound appears in 
stationary results [221 He] ). For reference on analytic results in the Newtonian case (p = 2) 
we refer to the introduction of [2]. 

Unfortunately, the Lipschitz truncation method of [T2| is not applicable to (ll.lOp since 
the system provides only control of dt{g'v) and not of dtv (unless p = const.) in suitable 
spaces. Alternatively we will use the L°°-truncation method, which was already successfully 
applied in [H EH EH] to prove existence of weak solutions for power-law type fluids if 


Throughout the paper we make the following assumptions: 

(Al) /: M ^ M is three-times differentiable such that there is some C > 0 satisfying 

|/"'(s)l < C(|s|+ 1) for all sGM 
and f”{s) > —a for all s G M and some a > 0. 

(A2) S: M X ^ satisfies 

|S(s,M)| < C(|M|P-i + 1) 

|S(si,M) - S(s 2,M)| < C|si - ssKlMr^ + 1) 

S(s,M) : M > a;|M|P - Cl 

for all M G s,si,S 2 G M, and some C,Ci,uj > 0, p G (^p^jCc). Here A : 

B = tr(A^B) and = {A G = A}. Moreover, we assume that 

S(c, •): ^ is strictly monotone for every c G M, i.e. Mi, M 2 G Rgy^ 

(S(c, Ml) - S(c, M 2 )) : (Ml - M 2 ) > 0 

with “ = “ if and only if Mi = M 2 . 
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(A3) q: M —7- M is twice continuously differentiable and strictly positive. Moreover, g, g', g" 
are bounded. 


We note that (Al) implies that 


\f{s)\<Ci\s\^ + l), 


\f{s)\<Ci\s\^ + l), 


\ns)\<Ci\s\^ + l) 


for all s G M and some C > 0. 

Our main result result is: 

Theorem 1. Let Q C be a bounded domain with -boundary, 0 < T < oo, and let 
(Al)-(A3) he valid. Then for every vq G (/?o G B 24 ^{LI) with dj^ipoloQ = 0, there is 

a weak solution 


V G L°°(0, T; Ll{n)) O LP{0, T; Wq tp G ^^’^(0, T; L‘^{n)) O L^(0, T; 


of (II.101) - (jl.Ill) in the sense that 

- / g{p)v ■ dtr}d{x,t) - / £»(v?o)vo • T 7 |t=o dx - / v {g{p)v + J) : Vr] d{x,t) 
J Q'j' Qt 

— R^-rid{x,t)+ S{ip,'Dv) :'Dr]d{x,t) = £o 0 Vy? : Vt 7 d(x, t) 


for all T] G C'°°([0,T) x 0)'^ with divT 7 = 0, (|l.lUI) -;- (|l.lUI) /i are satisfied pointwise almost 
everywhere and (ll.llD i . (jl.lip o. (jl.lip /i hold true in the sense of traces. 


The structure of the article is as follows: In Section [3] we prove existence of weak 
solutions of a suitable approximation of (|1.10l) - (ll.lip . where the convection terms pv • Vv 
and V • Vp are mollified. This is done with the aid of a (partial) Galerkin approximation 
for the Navier-Stokes part using unique solvability of (|1.10D q- (|1.10P /i for given v. Then in 
Section 0] we prove Theorem [T] with the aid of the L°°-truncation method by passing to the 
limit in the approximated system. 


2 Notation and Preliminaries 

For a, 6 G we denote a (8> 6 = ab'^ = G The standard Lebesgue spaces 

(with respect to the Lebesgue measure) are denoted by L^^M), 1 < p < oo, for some mea¬ 
surable M C LP{M;X) the A-valued Bochner space, and L^(0, T; A) = L^((0, T); A). 
The standard LP-Sobolev space is denoted by W'^'P{Ll). W^’^{Ll) is the closure of C^{Q) 
in and = kF™'’^(n), iL™(n) = W^’‘^{Ll). Finally, spaces of divergence free 

test functions will be denoted with a subscript “cr”. In particular, C^^(Q) = {tp G : 

div(p = 0}. L‘^{Q) is the closure of in Finally, Lq{Q) denotes the subspace 

of of functions wit zero mean value. 

By Bp q(]R‘^), s G M, 1 < p, g < oo, we denote the standard Besov space and its 

restriction to fl. We note that all results on interpolation of Sobolev- and Besov spaces for 
carry over to the spaces on Ll using a suitable extension operator, cf. e.g. the discussion 
in [U Section 2]. In particular, we have for every /c G N, 1 < po,Pi < oo that 

{w’^’P°{n),w’^+^’P^{n))e,p = B^+\n) - = ^ + -,fc g No, 

P Po Pi 
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for all 0 e (0,1), cf. [32l Section 2.4.2 Theorem 1]. Here (., denotes the real interpola¬ 
tion functor. 

Finally, the space of bounded and uniformly continuous functions /: [0, T] X for 
some Banach space X is denoted by BUC{\}}^T]]X). The space of weakly continuous 
functions /: [0,T] —)■ X is denoted by C'u;([0,T];X). 


3 The Approximated System 


In this section we consider the approximate system 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


gdt-y + (^'(Tgv) + J)) • Vv -I- - div S(</?, Dv) -|- Vtt 

divv 

dtip + (^ev) • Vip 


—Tg div (V(/j (8i V(/9), 
0 , 

A/i, 

- Aif 


together with 


(3.17) 


V = 0 
dj^(p = = 0 

v(0,-) = vo v?(0,-) 


on dQ X (0, T), 
on (912 X (0, r), 
(fio in 12, 


where R = — • V/i, J = —and Tg = e ^^Pa, where A = —PaA is the Stokes 

operator. Then £ 77^(12)*^ n 77q(12)‘^ n T^(12) for all v £ T^(12) and e > 0. Moreover, 
'kgf —>£^0 Pa^ strongly in L^(12)'^ for all f £ L^(12)‘^. We note that (|3.15p implies 


dtQ + div(f3'I'£V -|- J) = 7?. 


Therefore (I3.13P is equivalent to 

(3.18) dtig'v) + div(v 0 (^^('I'gv) -|- J)) — R^ — div S((/?, Dv) + Vtt = — Tg div (V(/7 (8) Vy?) 
Moreover, using 

- div(Vv; ® V^) = /rV(^ - V ^ 

we can replace — div(V(/9 (8> V(/9) on the right hand side of p3.13p by ^jX/(p, which will be 
used in the following. 

In order to show the existence of a solution of (I3.13|) - p3.16|) we first solve the Navier- 
Stokes and Cahn-Hilliard part separately. For this we need Besov spaces. However, they 
will not appear anymore in the remainder of the paper. 

Theorem 2. Let T > 0 and v £ L°°(0, T; 7y^(12)) he given. Then for every cpo £ 771 ^ 4 ( 12 ) 
with djg-gDo = 0 there is a unique ip £ 1F^’^(0, T; L^(f2)) n 7y^(0, T; 1F^’^(12)), which solves 
p3.15p - p3.16p together with p3.17|) 9 and p\t=o = ipo- Moreover, the mapping 

L°°(0, T; Ll{a)) 9 V ^ 5[v] = 99 £ 1F^’^(0, T; n L^(0, T; 

is continuous and hounded. 
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Remark 3. We note that the eondition tpo € is neeessary for the regularity of (f 

in the theorem because of 


cf. Amann 0 Chapter III, Theorem 4.10.2]. Moreover, because of B 2 ^{Tl) (7^(0), 

we have ip G L°°(0, T; (7^(0)). 

For the following we define 

= |rt G : J udx = 0^ for 1 < p < oo, 

= H\n) n 

and An - -H'(o)(f^) by 

— {ANipjtjj) jj-i jji = J Vp-Vifdx for all p,%lj & H^q-^{TI). 

As usual we identify u G with if i->- J^uifdx G Moreover, let Pog = 

5 — f^gdx. Then Pq: Lf‘{H) —)• L‘^{Q) is the orthogonal projection onto L^Q^(n). 

For the proof of Theorem [2] we will use: 

Lemma 4. Let T > 0 and v G L°°(0,T;L^(0)) be given. Then for every po S H^{Q) 
there is a unique solution 

p G L°°(0, T; H^{n)) n L^fd, T; H^{n)) n VFb2(0, T; R("J(S7)) 
of (j3.15l) - (j3.16l) together with (I3.17P 9 and p\t=o = Po the sense that 
dtp + '^e^ ■ Vy? = A^poit^x) in P'(0, T; (O)), 
where djpp\QQ = 0 and 

Po = -Ap + Pof'ip). 

Proof. In the following we assume without loss of generality that Jq po dx = 0. Otherwise 
we replace p and f by p — m and /(• + m) for m = Po(x) dx. 

The lemma will be a consequence of [H Theorem 4] and elliptic regularity theory. We 
refer to the Appendix |B] for a summary of this result and basic facts from the theory of 
monotone operators. To this end we define Bq = hI]Q^(fl), Hi = and A: P{A) —)> 

Ho, B: [0,r] xHi^Ho by 

A{p) = -An{-Ap +f]{p)), peV{A), 

V{A) = {p^ H^Q) n Rj,)(0) : d^p\an = 0}, 

B{t){p) = —'I/i;^{t) ■ Vp — aA^p for all G Ri, t G [0,T], 
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where f{s) = /o(s) — with /o: M ^ M convex and a > 0. Moreover, let Eq : —)■ 

(—oo,+oo] be defined by Eo{if) = ^ /o(¥^)dx if ip £ dom(£'o) := and 

Eo{p>) = +00 else. Because of Corollary [TT] in the appendix, A = Now [TJ 

Theorem 4], cf. Theorem [9] in the appendix, implies the existence of a unique solution 
ip € L°°{0,T; n W^’'^{0,T; with ip{t) G E>{A) for almost every t G (0,T) 

since 

\\^{'t)\\c{Hi,Ho) ^ \\'^e^\\L’=°{ 0 ,T-,H^) + « < C'e||v||Loo(o,T;L 2 ) +« 

for all most every t G [0,T]. Using elliptic regularity theory and the equation, one obtains 
additionally </? G L^(0, T;□ 

Proof of Theorem [M' Due to Lemma 0] it remains to prove higher regularity of the solution 
(p . To this end we use that 

dtip + AjvAip = Po(f"((p)Aip + f"'((p)lV(pl^) - TgV • Vcp =: E, 
in L^(0,T; where 

||/"((^)A(^||i2(gj,) < ||/"((^)||x,oo(o,T;L3)II^V^IIl2(0,T;L6) < OO 

ll/"'(‘^)|V^niz,2(2^) < ||/"'(^)||L-(0,T;La)||V99||i4(o,r;L6) 

^ C{\\ip\\L°°{0,T;Hi) + 1)II‘^IIl4(o,T;//2) < OO 

llTeV • Vy5||i2(Qj,) < ||'I'eV||x,oo(gj,)||Vy5||i2(Qj,) < Ce || v|| ^oo (o,T;L 2 ) || (/3|| j ;^2 < OO 

since |/"(s)| < ^(s^+l), |/'"(s)| < C(|s|+1) and yp G L°°{0,T; H\n))nL‘^{0,T; H^{n)) ^ 
L‘^(0, T; H‘^{Q)). Hence, by linear theory we obtain ip G L^(0, T ; T; L^(D)) 

due to m Theorem 2.3]. In particular ip G L°°{0,T; Using similar estimates as 

above we obtain E G L^(0, T; L^(fl)). Therefore applying |15l Theorem 2.3] once more, we 
obtain ip G T; L‘^{n)) n L^(0, T; H^{n)). □ 

Lemma 5. Let k > 0 be arbitrary. Then the solution operator iS[v] from Theorem 
satisfies 

sup ||5[vi] -5 [v2]|Lo < c{k,S,T)T^ sup ||vi - V2|Lg 
(o,r) ’ (o,T) 

for any vi, V 2 belonging to the set 

Mk = {u G L°°(0,r;L^(D)) : ||u||£,oo(o,r;L2(n)) < 

Here c{k, T) is non-decrea.sing with respect to T > 0. 

Proof. Let vi, V 2 G be given. We set ipi = 5[vi] and ip 2 = 5[v2] such that Theorem [2] 
implies 

(3.19) II¥^i||l°°(Qt) + II¥^2||l°°(Qj.) < c{k). 

First of all we have 

dfiipi - ip2) + (^eVi) • Vipi - (TeV2) • V ip2 = A{f'{ipi) - f'{ip2)) “ A^ {ipi - (^ 2 )- 
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Testing with ^pi — (p2 we obtain 


^-4 [ Wi-^2\^dx+ f |A((/?i - (^2)Pdx 

t/ J O 


= - j (^(TeVi) • V(/?i - (T£V2) • V^p2^ {(pi - (P2) dx 

+ [ - f{^2)) ■ A((/?1 - (P2) dx =: Ji + J2. 

Jn 


The first term we estimate via 


Jl = y (^(TeVi) ipi - (T£V 2 ) ip 2 ^ • V(<y 9 i - (^2) dx 


< 


(TeVi) (^1 - (T£V2) (/?2 dx+ \V{ifi - ip2)\^ dx 

Jn 

dx + c{ 6 ) I Iv^i - (p2\^ dx + 5 I I A(<y 9 i - (^2)^ dx 
Jfi Jq 


< / (^£Vi) (fi - (T£V2) (f2 

JQ 

= : Ji + J 2 +j 3 ^ 

where (5 > 0 is arbitrary. On account of vi E (| 3 . 19 l) and the continuity of Tg we have 

< 2 / e^i) {ipi - (p2)f dx + 2 f |Te(v2 - Vi)<y92|^dx 


< 2 ||^'£;Vi||^ / \ipi - ip2f dx + 2\\ip2\\oo [ |T£(v 2 - Vi)|^dx 

J Q Jq 

< c{K,e) / \ipi - (p 2 ^ dx + c{k) / |v 2 -vi|^dx. 

Jo, Jo, 

Moreover, we obtain by (A3) and Young’s inequality 

J2= / f”{(p2 + 9{(pi-‘P2))de{ipi-(p2)^{‘Pl-‘P2)dx 

JQ Jo 

< c{k) / |(^1 - (/ 92 ||A((^i - ^ 52)1 dx 

Jn 

< c{6,k) / |(/?i - 992 !^ dx + 5 / |A( 99 i - 9 ? 2 )|^ dx 

Jq Jq 

Plugging everything together we have shown 

< c{5,k,s) / I 991 - 992 !^ dx + 2(5 / |A(99i - 992 )!^ dx + c(k) / |v 2 -vi|^dx, 

Jo Jo Jo 

for any (5 > 0. Now we choose (5 small enough and integrate with respect to t such that 
(note that ipi and 992 have the same initial datum) 

[ I 991 - 992|^dx+ [ [ |A(99 i - 992 ) 1 ^ d(x,t) 

Jo Jo Jo 
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< c{K,e) 



The claim follows by Gronwall’s lemma. □ 

Theorem 6. Let Assumptions (A1)-(A4) hold true and let e > 0, vq G Then there 

is a unique solution (v, ip, p) to (j3.13p - (j3.16p and (j3.17p i - (j3.17p -i with v G L°°{0, T ; L^(n))n 
L^’(0, T; ip G T; L‘^m n T; W^^\Ll)), // G 7.2(0, T; 

Proof. The aim is to find a function v G L°°(0,T; (Q)) n U’{0,T-,Wq’^{LI)) such that 


(3.20) / Qv ■ dtrid{x,t) + / S((/9, Dv) : Dt 7 d(3:, t) 

j Qj^ j Qrp 


/ ^v(8)'I'£(v) : VT 7 d(x,t) + / df ^ div{Vip'S) V(p) ■ rid{x,t) 

J Qj- J Qj- 

i J : VT7d(x,t) + / £»o7^o • ^(0, Odx 
Jn 


' Qt 

+ f ■ r]d{x,t) + f 

J Qrp J Qrp 


for all T] G C'^([0, r) x Q), where we abbreviated J = —R = —• V/x and 
g = q{(p) and ip is the solution of ()3.17l) ':i- (l3.17P /i due to Theorem [2] and p is defined by 
(j3.17p .i. Note that we used (j3.18p . 


Step 1: finite dimensional approximation. 

We separate space and time and approximate the corresponding Sobolev space by a finite 
dimensional subspace which leads to a system of ODEs (Galerkin Ansatz). To solve this 
we follow the approach for compressible Navier-Stokes equations introduced in [T8l Section 
2]. From |29l Appendix] we infer the existence of a sequence (A^) C M and a sequence of 
functions (w^) C V17 q’^(12), I G N, such that 

i) Wfc is an eigenvector to the eigenvalue in the sense that: 




Afc / Wfc • T7 dx 


for all r] G 


ii) WfcWm dx = Skm for all k,m G N, 
hi) 1 < Ai < A 2 < ... and A^ —)• 00 , 

^ A:,m G N, 

I 2 

v) (wfc)fcgH is a basis of 1 Tq]^(0). 

We choose / > | + 1 such that Wq‘^{LI) ^ lT^’°°(n) We are looking for an approximated 
solution of the form 


N 

E N 

Ck Wfc 

k=l 
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where : (0,T) —)■ So we have for a.e. t e (0, T) 


VAr(f) £ Xn := span{wi, ...,wn}. 


We will construct such that is a solution to (here we denote go = g{(po)) 

/ Q{(PN{t))vN{t) ■'Wkdx = / £»oVo-Wfcdx+ / VAT (g) ^)((^Ar)(^£VAr + J) : Vwfcd(x,CT) 
</ J ^ J Qt 

(3.21) +[ Rjsf'^-Wkdix^a)- [ S((/9Ar, Dvat) : Dw^ d(x, ct) 

jQt ^ JQt 

— / 'he div (V</?Ar (g) V(/97 v) • Wfc d(x, (t) for all t£[0,T], 

JQt 

for k = 1, where 


(^Ar=5[vAr], Jat = -^(</?Ar)V/iAr, 


dg 

Rn = 


with fj.N = f'{^N) — A(/77v. Here S denotes the solution operator from Theorem [2l We 
introduce the operator M.i\j{q\ '■ Xm ^ X'j^ with 


(Wl7v[^>]u,w) 


• w dx 


for all 


u, w £ Wat. 


As g{i^) is bounded from below by some g > 0 the operator Mm\q] is invertible and we 
have 


(3.22) \\Mn[ 0^] ^-MNie"^] ^ 

So equation (|3.21l) is equivalent to 

(3.23) VAr(t) = AdAr[£>(5[vAr])] ^ (^{go^oY + J ^[^[vAr], vat ] da^ , 

Here (^qVo)^ denotes the unique A^y-representative of £>oVo and we abbreviated 


{M[ip,VN],w) = 


J = 


/ Vat (g {g{ip){df^VN) + J)) : Vwdx + / R^ ■ wdx 
Jn Jn ^ 

/ S((/;, Dvat) : Dw d(x, (t) — / 'I'g div (Vy? (g) V(/?) • w dx, where 

</ £1 ./si 

-|^(^)V/r, i? = -v|^((^)-Vp, gt = r{ip)-Xtp. 


As the space X^ is of hnite dimension it is not hard to see that a local solution vat £ 
C'([0,Tat]; Aat) to (I3.23P can be found by a standard fixed point argument provided Tat 
is small enough (recall (I3.22p and Lemma [5]). Then V[^(pn'v^n) £ C*^([0, Tat]; Aat) as a 
consequence of (I3.2ip . Here denotes the orthogonal projection Wq^(H) —)• Aat. In 
order to gain a global solution on [0, T] we have to show that vat stays bounded in Aat. We 
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differentiate (j3.21l) in time and use v^r as a test-function. Moreover, we use ipN = 5[vAr] 
as a test-function in (j3.15p and dt^pN in (|3.16p . Then we obtain by (A2) 


sup f [ ^i7v|vArP+ / dx+ [ f{pN)dx] 

o<t<Tjv \Jn Jn ^ JQ J 

rTff !■ pTff I- 

+ / |Dvw|^dxdt-F / / dxdt < C(yo,po), 

Jo Jn Jo Jn 


where CPYg, </?o) does neither depend on T/v nor on N. Hence we gain a global solution in 
time, i.e. = T and uniform a priori estimates in N. In fact, as g is strictly positive 
(uniformly in N) we have 


(3.25) 


sup / |v^(t,-)l^dx+ / |Vv^|Pd(x,t) < C 
tG{0,T') J Q. J Qj' 


for some C > 0, where we also used Korn’s inequality. By Theorem [2] and Remark [3] there 
holds 


(3.26) PN G T^(0,r;iy^’2(0))mT^’^(0,r;L2(o)) L°°{0,T-,C\n)) 

as well as 

(3.27) pN e L\0,T-,W^’\n)), 

both uniformly in N. For the last statement we used equation (13.151) together with (I3.26P 
and TgVAr • VipN G L°^{Qt) (combining p3.25p and Tg: T^(H) —^ ^ 


Step 3: weak convergence. 

On Passing to a subsequence we gain from (I3.25p - (|3.27l) 


(3.28) 

VAT - 

V 

in 


(3.29) 

VAT - 

V 

in 

L^{0,T;L^{n)), 

(3.30) 

Pn - 

p 

in 


(3.31) 

Pn - 


in 

iyi’4(0,r;L2(o)), 

(3.32) 

Pn - 

p 

in 



as —>• oo. The convergences of ip m above imply in particular 


(3.33) 

VifN V 99 

in L^{Qt), 

(3.34) 

ifN ^ ill 

L'^( 0 ,r; 0 ^( 11 )) 

(3.35) 

Jjv — y J in 

L\Qt), 

(3.36) 

R]\j —)* R in 

L\Qt) 


as —?> oo, where g < oo is arbitrary. Hence we can pass to the limit in the Cahn-Hilliard 

equation and obtain 


dtp + ('I'ev) • V(/9 = A/i a.e. in Qt 
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for /i = f'i^p) — A(p together with 


dj\fip = = 0 

^9(0,-) = ipo 


on X (0, T), 
in 0. 


In order to pass to the limit in the convective term we need compactness of vjv in LP‘{Qt)- 
We obtain from (j3.2ip 


' Qt 


' Qt 

+ / 'h, 
J Qy 


dtVj^iQN^N) • J? dx = [ dtiQN^N) ■ V^^ri dx 
J Qj, 

f S{(pN, Dvat) : DVi^r] d(a;, t) + [ Qn'vn ® ^e{^N) ■ d(x, t) 

J 0/p J Qrp 

■'P[^rid{x,t) + [ Rn^ ■Vi^rid{x,t) + f vat 0 Jat : VP/^T7d(x,t) 

J Qrj^ J Qj^ 

= -[ S{(Pn,'Dvn) ■'^'P[^r]d{x,t) + [ qn^n ^ei^N) ■d{x,t) 

J Qrp J Qrp 

- j VA-^^^,{^ 1 NV^N) : ViPfr?d(x,t) - f VA^\Rn^) : ViP/^r?d(x, t) 

J Qrp J Qrp 

+ [ VTV® Jjv : VP;%d(x,t) =: / Hjv : d(x, t) 

J Qrp J Qp 

for all 77 G L^(0, T; Wq ^(n)). Here denotes the orthogonal projection into Xn with 
respect to the Wq’ (H) inner product. Note that due to the choice of wi, ...,'wn the pro¬ 
jection is orthogonal with respect to the L^(H) inner product as well. Moreover, the 
operator is the solution operator to the Laplace equation on H with respect to zero 
Dirichlet boundary values. On account of (|3.25l) - (l3.4ip we have uniformly in N 


(3.37) 

(3.38) 

(3.39) 


gTVVTv ® ^'s(v7v) G 7.°°(0,T; l2(H)), 
^'^(/iTvVy.Tv) GL^( 0 ,r;L°°(H)), 
i?,v,JivGL^( 0 ,r;L''(‘^)(L!)), 


where q{d) = if d > 3 and q{d) < 00 arbitrary if d = 2. Here we also used the properties 
of d/g. So we have uniformly in N 

(3.40) Hjv G Lp'{Qt) 

using the properties of and p > 

Step 3: strong convergence. 
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We get on account of (I3.40p and Sobolev’s embedding (recall the choice of 1) 
Wdt'Pl {QN'Vn)\\ip> — \\dtVl 


(3.41) 


< 


< c 


sup 

[ dtiQN^N) -Pi^^dix, 

.t) 

^^‘^^Lp{0,T;WQ’^(n))<l 

Qt 


sup 

[ Hat : VP;'^</5d(x,t) 


ll<p|| 12 J 

"^"rp(o,T;wy^(n))<i 

1 

Qt 

^ f / \P / 

r 

sup 

/ |H^rd(x,t) 


ll<pll 12 ' 

"^"rp{o,T;wy^(n))<i 

V JQt / V 

J Qrp 


\VPi^<f\Pd{x,t) 


sup 


V^\\LP'{0,T-,L°°{n)) — ^■ 


Combining (13.2811 . (I3.34|l and (I3.4ip with the Aubin-Lions compactness Theorem shows 
Pf^( qn'^n) —^ PaiQ"^) in T^(0, T; L^(14)). Note that we have the pointwise convergence 
—)• P(j in L^{Q). Due to this and (I3.29h we have 

I |v^VArpd(x,t) = / Vj^iQN^N) ■^Nd{x,t) ^ / Pa{Q-v) ■vd{x,t) 

J Qrp J Qrp J Qrp 


' Qt 


|^vpd(x,t) 


for N —)■ oo. This means that -y/^vjy converges strongly in LP‘{Qt). As qn is strictly 
positive we have compactness of v^r as well and hence for N —)• oo 


(3.42) 

(3.43) 


vtv —^ V in L^{Qt) for all s < p 


d + 2 
d ’ 


vtv —)• V in L'^( 0 , T; L^(D)) for all q < oo. 

,d +2 


Here we have used that L°°(0, T; ^^(D)) n LP{0, T; Wpi(D)) LP—{Qt). This and (iTMIl 
yield together with the continuity of 

(3.44) ^/^VN in L^{Qt), 

(3.45) ^tvvat (8)'I'eVTv —)• (g) in L^{Qt), 

for all s < Due to (I3.25P and (A2) we know that S(v9jv, D(vAr)) is bounded in L^'(Q) 

thus 

(3.46) S((/?Ar,D(vAr)) ^ S in V'{Q). 

Finally we can pass to the limit in (I3.21|l such that 

/ Q{p{t))^{t) ■ rjdx = / ^)oVo-T 7 dx+ / / v (g) ^('I'gV + J) : V 77 dx du 

JQ. Jn Jo Jn 


(3.47) 


+ 



■ rjdxda — 


/o Jn 

rt 



S : D 77 dx da 


0 Jn 


n 'kg div (V (/9 (g) V(/j) ■ rjdxda 
! 
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for all T} G C'^(n) and t G (0,r). This implies the following equation for dt{Q^) (using 
continuity of Tg) 


' Qt 


■ dt'nd{x,t) = / £ioVo • ^7(0) dx + / w ® + i) ■.Vr}d{x,t) 

J Q. J Q’j’ 

(3.48) + f ■ r]d{x,t) - j S:DT 7 d(x,t) 

J Qrp J Qrp 

— / Tg div (V (/9 ( 8 ) Vy?) • 77 d(x, t) 

J Qrp 

valid for rj G C^([0,T) X £7). By density argument this yields 

dt{g^r)eLP\ 0 ,T-,W-^’P'm. 

Step 4: monotone operator theory. 

We apply monotone operator theory to show 

(3.49) S = S((/?,Dv). 

To this end we have to study the term 

j (s{(pN, Dvat) - S((/?, Dv)^ : D(v 7 v - v) d(x, a) 

= [ (S - S((/?Ar,DvAr)) : Dvd(x,o-) - S((/?, Dv) : D(vAr - v) d(x, u) 

jQt ^ ^ JQt 


+ 


/ S((/9W) Dvjv) : Dvat d(x, t) — / S:Dvd(x,(T). 
Qt Qt 


We get from (j3.28jl and (13.461) 


/ ( S — S((^7V) Dvat) ) : Dv d(x, a) —> 0, N ^ 00 , 
JQt 


'Qt 


S{(p, Dv) : D(vAr — v) d(x, a) —> 0, N ^ 00 


for all t G (0,T). Moreover, by (I3.2ip (differentiated with respect to t) and (I3.48P there 
holds 

/ S((/J 7 v, Dvv) : Dvtv d(x, cr) — / S:Dvd(x,fj) 

JQt JQt 

= \ j ^’iv(0)|vAr(0)pdx-i f £io|vopdx+ f • 77 d(x, t) 

J 0 J ^^t 

— / Te div (V(/?Ar ® V(^Ar) • vw d(x, u) + / div (Vy? 0 V<y9) • vd(x, u) 

JQt JQt 

[ QN{t)WN{t)\^dx+^ f £<(t)|v(t)|2dx 

0 J 

=: (O)+ (/) + (//)+ (///). 
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The first term can be treated using 

/ t'Ar( 0 )|vAr( 0 )|^ dx = [ V^iQovo) ■ MN[Qo]~^i'Pi^{eo^oy)dx —;► [ ^lolvol^dx, N^oo. 

J ^ J ^ ■/ o 

Here we used strong convergence of P/^(£»oVo) and pointwise convergence of the operator 
-Ad^M^'o]- So 

lim (O) = 0. 

N^oo 

On account of (|3.36l) and (I3.42p there holds 

lim (/) = 0 . 

N^OO 

We deduce from (I3.28p . ()3.34p and continuity of that 

lim (II) = 0. 

TV^oo 

Finally, we obtain by (I3.44p for a.e. t G (0, T) 

lim (///) =lim f gN{t)\^rN{t)\'^ dx + ]- [ g{t)\^r{t)\‘^ dx = 0 

N^oo 2 N^oo 2 Jq 

at least after taking a subsequence. Hence we end up with 

(s{ipN, Dvat) - S(v 3 , Dv)^ : D(vAr - v) d(x, t) —0, N ^ 

for almost every t G (0,T). We have by assumption (A2) for any 1} G (^, l) 

{S{ipN, Dvjv) - S((/?, Dvtv)) : D(vAr - v) d(x, a) 

C [ |(/?Ar - + |DvAr|P"^)'^(|DvAr| + |Dv|)’^d(x,cr) 

Jo.t 


oo 


’Qt 


'Qt 


< 


< C { I \(PN - ‘fl 1-’’ d(x,(j) 


1 - 1 ? 


'Qt 


(1 + |DvAr|^ + |Dv|P) d(x, a] 


for N ^ oo, where we used ()3.34p and (I3.25p . This finally implies using monotonicity of S 

J (^(S((/?,DvAr) - S(¥5 ,Dv)) : D(vAr - v)^ d(x,cj) —^ 0, N ^ oo 

for almost every t G (0,r). By monotonicity of S (which follows from (A2)) we obtain 

(1091) . □ 
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4 Proof of the Main Result 


In the following let (v„, (/?„, //„) be a solution to (I 3 . 13 p - (l 3 . 17 p with e = ^ which exists due 
to Theorem [6l So there holds 


( 4 . 50 ) - 


/ Qn^n-dt-nd{x,t)+ S(<y9„,Dv„) : DT7d(x,t) 

J Qrp J Qrp 

= (g) Ti/„(v„) : Dr7d(x,t) + / • ry d(x, t) 

JQt JQt 


+ / Rn^ ■ 'nd{x,t) + / (g) J„ ; D?7d(x,t) + / •’iCO, •) dx 

J Qjp J Qj- J Q 


for all r] G C^^([ 0 ,r) x 0 ). Here we have 3n = -|^V/Xn, = -V|^(y?n) • V/Xn and 

— Q{}Pn)- 


Step 1: weak convergence. 

We obtain the following a priori estimates (testing the momentum equation with v„, the 
equation for tpn with and the equation for with dtipn) 

sup ( [ ^„|v„pdx+ [ ■^^^^dx+ [ f{ipn)dx'] 
te(o,T) \Jn Jn ^ Jn J 

+ / |Vv„|Pd(x,t) + f |V/r„|^d(x,t) < C. 

J Qj^ J Qrp 


This implies the following convergences (after choosing appropriate subsequences) 

v„-v in LP( 0 ,T;Wi’^(O)), 

in L°°( 0 ,r;L 2 (H)), 

S((/ 7 n,Dvn) ^ S in Lp'{Qt), 

(4-51) _ ^ 

Vn < 8 ) Qn^l/ni^n) ^ H m 2 d (Qj,) 

(fn^^ in L^( 0 ,r;C'^(H)), 

Vifn Vif in L^(Qt) 

as n —>• oo. By Theorem[ 2 ]and Remark[ 3 ]we gain due to the boundedness v„ G L°°( 0 , T; L^(H)) 
ipn G L^( 0 , T; p T; L^{n)) ^ L°°( 0 , T; C\n)) 

uniformly in n and therefore 

\\^n\\L^{0,T]W^'^) + \\f{^n)\\L^{0,T]L^) < C. 


So, hnally we have 


( 4 . 52 ) 


in L°°(Qt), 
in L<1{0,T-C\II)), 
in L\0,T-W^’\n)), 
Q{‘Pn) q{^) in L°°{Qt) 
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as n —)■ oo for all q < oo. As we have dt^fn S L‘^{Qt) and v„ • V^pn S L'^{Qt) uniformly in 
n (recall (I4.5ip i and ()4.52p 9l. equation (I3.15P yields 

Apn G L^{Qt) 

uniformly in n. Combining this with the a priori estimates we also have 

(4.53) e L\Qt) 
uniformly in n and hence by (A3) 

(4.54) VJ„ e L\Qt) 
uniformly in re. From (I4.50p and the a priori estimates we gain 

(4.55) dtPaiQn^n) G L"°(0,T;fF-i’^“(O)) 

uniformly in re for some sq > 1 (not that we have only control over solenoidal test-functions 
but div(£»„v„) / 0). In addition, we have by (|4.52p i and the continuity of Pa on VF^’^’(n) 

(4.56) VPaiQn^n) G ^^(Qt). 

Both together yield compactness of PaiOn'^n) in IP{Qt) by Aubin-Lions compactness the¬ 
orem. So we have 


(4.57) PaiQn^n) -^n^oo Pu{Q^) in L^{Qt) 

using (I4.52p . Due to (I4.57P and (|4.5ip9 we have 



lv^v„l^d(x,t) 



Pai,Qn^n) 



V d(x, t) 



|^vpd(a:,t) 


for re —>• oo. This implies that converges strongly in L?‘{Qt). As is strictly 

positive we have compactness of as well and hence 


(4.58) Vn ^ V in L^{Qt) for all s , 

a 

(4.59) v„ —>■ V in L'^(0, T; L^(D)) for all q < oo. 


Moreover, we have 
(4.60) 




J = -gV/r in L\Qt), 


Jn ^ 

Vn <8) Jn ^ V (g) J in L'‘°{Qt). 
for some sq > 1. On account of (I4.54p we have 
(4.61) V(v„ (g) J„) ^ V(v (g) J) in L^°{Qt) 
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combining (I4.54p with ()4.58p . 

Combining ()4.51l) /i and (j4.58l) yields H = v (8) gv. This means we also have the conver¬ 
gences 


d+2 


(4.62) 


V(v„ (g) £»n'I'i/n(vn)) ^ V(v(8) £»v) in {Qt), 


using the pointwise convergence '^i/n —^ Pa- We have the limit equation 

/ Q'v ■ dtri(l{x,t) + / S:V77d(x,t)— / v (g) : Vry d(x, t) 

j Q rj. j Q rp j Q rp 

V (g) J 


(4.63) 


St d Qt 

' R^-rid{x,t)+ / V (g) J : Vt 7 d(x, t)-|- / fiVip ■ rid{x,t) 
Qt Qt Qt 


for all T] S C^^{Qt) and the equation for the difference Qn'Vn — Q'v 

- {QnVn-gv)-dt'nd{x,t)+ {S{^Pn,'DVn) - s) : Vr]d{x,t) 
J Qrp J Qj’ 

- (v„ (g) £>„Ti/„(v„) - V (g) £iv) : V77d(x,t) 

JO.T 


(4.64) 


'St 


[Rn^ - R^) ■ V d(a:, t)+ / (v„ (g) J„ - V (g) J) : V 77 d(x, t) 


'Qt 


+ (^l/n(/^nV</?n) - /iVyj) • r/d(3:,t). 

J Qrp 


Step 2: pressure function and strong convergence. 

Now we define 


H”:=H- + H-, H^=S(^„,Dv„)-S, 

H 2 := -V„ (g) Qn^l/n{^n) + V (g) £»V - (v„ (g) J„ - V (g) j) 

where is the solution operator to the Laplace equation on 14 with respect to zero 
boundary conditions. Then we have 

-[ {Qn'Vn - g'v) ■ dtrid{x,t) = - f Ui :V-nd{x,t) + [ div ■ rj d{x,t). 

J Qrp J Qrp J Qrjp 

Moreover, we know 

(4.65) in Lp'{Qt), 

(4.66) H 2 ,divH 2^0 in L^°{Qt) for some sq > 1 

and introduce the pressure in accordance to Theorem [7] and Corollary [8l We obtain func¬ 
tions tt)), tt”, 7r2 such that for some r G (1,2] and any 14'd 14 

lkhllL°°(0,T;L’-(O)) < C n|H"||Lr(Qy) -h IIu"-11^00(0,T;L2 (o))] , 
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IIlp'(Qt) - IIlp'(Qt) 


< cIlH’i 


F2 IIl’'(St) - 


< cIlH 


2 IIL’-(Qt) 


ll^^2 llL’'((0,T)xf2') < c(r2') IIVH 2 llLr-(gy) 
for some c > 0 and 

- {gnVn-Q^)-dtr]d{x,t)+ 7r^dtdiyr]d{x,t) 

JOrr JO.'r 


(4.67) 


/ (H” - 7r"l) : V77d(x,f) + / div (H 2 - tt^/) • ?? d(x, f) 

J Qrp J Qrp 



0 

in 

L°°(0,r;L’'(S2)) 

^1 - 

0 

in 


vr2"- 

0 

in 

L^{Qt), 

Vtt^- 

0 

in 

L^((0,r) X !2') 


for all T] S C^{Qt)^- We have the following convergences for the pressure functions: 


(4.68) 


A main difference to the approach in [33] is that tt)) is not harmonic. In fact, we have by 
Theorem [3 a) 

(4.69) Avr)) = -div(£inV„- qw) = -V£>„ • v„ + • v G L°°(0, T; L^(n)) 

uniformly in n (recall (I4.5ip and ()4.52p l. Moreover, (|4.52p and ()4.58p imply 

d + 2 


(4.70) 




0 in L^{Qt) for all s < p- 


d 


For any Q' d 12" d 12 and any q,r G (1, 00 ) we have 


(4.71) 


IQ' 


dx < c{q,n',n”)( / | A7r;)(t)|'^ dx + ( / K(t)|'’dx 


IQ" 


JQ" 


This estimate is a consequence of |23l Theorem 9.11] in combination with Sobolev-embeddings. 
We gain on account of (|4.68I) - (I4.70I) and (I4.7ip 


(4.72) for all s<p^, 

uniformly in n. Now we are concerned with compactness of vr)). We set := Qn^n — Q'v + 
Vtt)) and have 

(4.73) q„GLP(0,r;W/,f(O)) 

uniformly in n (combining (|4.5ip . (|4.52p and (|4.72p ). Moreover, there holds by (|4.67p 

(4.74) q„G W^’''(0,r;W-i’^(I2)) 
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uniformly in n. This implies by the Aubin-Lions Theorem using (I4.72p 


(4.75) 


0 in L'^(0,T;L? (n)) for all q<p 


d + 2 

IT' 


Combining this with ()4.58l) and (j4.52P /i yields 


(4.76) 


V7r)( —)• 0 in L'^(0,T;L^^^(n)) for all q <p 


d + 2 

d ■ 


We can assume that '^h(x, t) dx = 0 for a.e. t, such that 


(4.77) 


0 in L'^(0,T;L^^^(n)) for all q<p 


d + 2 


In order to show (j4.77p we make use of the Bogovskh operator introduced in [TO]. It is a 
solution operator to divergence equation with respect to zero boundary conditions. Taking 
an arbitrary ball B Q and the Bogovskh operator Bog^ with respect to this ball we have 
Bogg : L++B) —)■ Wq'^{B) for all r € (l,oo). So we gain 


TT 


hWL+B) 


V&C^iB) 


V&Cg°{B) 


< sup 
V&Cg-iB), II 

< sup 
V&Cg-(B), II 


sup / T^hid- Vb) dx 

B), \\v\\^<=lJn 

sup / 7r)( divBog 5 (r/ - r/s) dx 

B), ||,7||„,=1 Jo 

/ V7r;)-BogB(??-??B)dx 


,=1 JB 

l|V7r;(||i,(B)|| Bogsiv - Vb)\\li'(b) 


„'=i 


< cW'^'XhWl+B)- 

Integrating in time and using ()4.76h yields (I4.77h . Finally (I4.7nh - (l4.77p imply 


(4.78) 


7 r)( —)■ 0 in L'^(0, T; for all q < p 


d + 2 


In the following we need to show that S = S(<^, Dv). 

Step 3: L°°-truncation and monotone operator theory. 

By density arguments we are allowed to test with r] G LP(0, T; IFQ^’^(n)'^) n L°°{QtY 
in (|4.67l) . Since the function v does not belong to this class, the L°°-truncation is an 
appropriate method (see [20| for the steady case and [33| for the unsteady problem). 

We define and Hl, L G Nq, by 

hL{s)-.= [\Li0)0de, Hl{^) := h^m), 

Jo 

L 

Tl:=^V’ 2-G V’5('S) := V'('5s), 

£=1 
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where il) € C°°([ 0 ,oo)) with 0 < -i/) < 1 , -0 = 1 on [ 0 , 1 ], -0 = 0 on [ 2 ,oo) and 0 < —' 0 ' < 2 . 
Now we use in ()4.67p the test-function rj = r/T 2 ,(|q„|)q„, where rj S C^{Q.). This yields 
for a.e. t G (0,T) using qn(0) = 0 (recall Theorem [7] a)) 

[ 7]HL{qn{t))dx = [ ((9tqn,?7TL(|q„|)q„) du 

Jn Jo 

= - [ f r/(H” - </) : V(TL(|q„|)qn) d(x,cr) 

Jn Jo 

-[ [ (H”-tt"/) : Vr/(g)TL(|q„|)q„d(3:,fj) 

Jn Jo 

+ [ [ 7/div(H^-7r^/)Ti(|q„|)q„d(x, a) =:(/) +(//) + (///) 

Jn Jo 

and dtcin G L^'{0,T-, n L°°(n))'), recall (I4.74p . The aim of the following observa¬ 

tions is to show that the terms {II) and (HI) vanish for n —)• oo which gives the same for 
(I) (note that the term on the left hand side vanishes for a.e. t by p4.75p l. By monotone 
operator theory this yields Dv„ —)• Dv a.e. Due to the construction of we obtain, after 
passing to a subsequence, 

(4.79) Tidq^Dqn— ^0 in L'^{Qt) as n ^ oo 

for all r < oo (first, we have boundedness in L^, then the strong convergence follows from 
(j4.75p l. This implies 


(//),(///) ^ 0, n^oo, 

as a consequence of (I4.65P and p4.66p . Plugging all together, we have shown 

limsup f 7?(S(¥5n,Dv„) - S)) : Ti,(|q„|)Dq„d(x,u) 
n^oo J 

(4.80) ^limsup/ (-r?) (S((/?n, Dv„) - S)) : VTi,(|q„|) (g) q„ d(x, a) 

n —>00 J 

-h limsup / T/TT^'div (Tx,(|q„|)q„) d(x,cr). 
n—>-oo J 

Now we want to show that the right-hand-side is bounded in L G N. Since divq„ = 0, 
there holds 

limsup / T/TTi div (TL(|q„|)q„) d(x,fj) 
n^co J 

= limsup/ T/7r”VTL(|q„|) • q„d(x,(T). 
n—>cxD J 

So, by (I4.68p . we only need to show 

(4.81) VTi(|qn|)qn G {0, T; Ll^{n)) 
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uniformly in L and n to conclude 


(4.82) limsup / 7/(S(<y9„, Dv„) - S)) : TL(|q„|)Dq„ d(x, a) < if. 

n—>cxD J Qf. 

We have for all i € No 

|v{V' 2 -Hlqnl)}qn| < ® 

< -2“dqn|^'(2“dqn|)|Vqn| < c\Vq^ri\XAt, where 
:= { 2 ^ < |qn| < 2^+^} . 

This implies 

L L 

|VTi,(|q„|)q„| < |V{V' 2 -Hlqn|)}qn| < |Vq„|xA^ < c|Vq”|. 


e=o 


£=0 


This yields (I4.8ip and hence also (I4.82p is shown. Now we consider 


SL,n : = 


[ r]{S{(pn,Dvn) - S)) : TL(|q„|)Dq„ d(3:, u). 

jQt 


On account of (j4.82p we have < K independent of L and n. Thus, using Cantor’s 
diagonalizing principle we obtain a subsequence with 

(^i,nk ■= / ??(S(</?nfc,Dv„J - S)) : V' 2 -H|qnJ)Dq„j^ d(3:,fj )—> ai as I 00 
JQt 

for all i G No- We know as a consequence of the monotonicity of S assumed in (A2), 
Vpn -A Vp in L^{Qt) for any s < 00 , and —)■ 0 in for n — 00 , see (I4.78h . 

that (Ti > 0 for all i G N. Moreover, a£ is increasing in £. This implies on account of (|4.82p 

„ . ao + ai + ... + ai K 

0 <<T„< -J- <J 

for all £ G N. Hence we have cjo = 0 and therefore 

/ r/(S((/?„,Dv„) - S)) : V’i(|qn|)Dq„d(x,t) —)■ 0, n-)■ 00 . 

J Qjt 

Due to (I4.5ip . (j4.52l) . (j4.58p uniform boundedness of 'i/'i(|qn|) and (14.781) we have for n —00 

[ p{S{ipn,'Dvn) - S((/?,Dv)) : V’i(|qn|)DV7r)(d(x,fj) —^ 0, 
jQ.t 


'Qt 


p{S{(pn,'D^/n) - S(vJ,Dv)) : V’l(|qn|)D((£'n “ Q)^n) d(x, u) -^ 0, 


and hence 


/ ? 7 (S((/?n,Dv„) - S(¥5,Dv)) : V'i(|qn|)D(^)(v„ - v))d(3:,fj) —^0 as noo. 

JQt 
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Using Vg G L°°{Qt) (14.5111 and (14.5811 this implies 

(4.83) [ OT(S((/?n,Dv„) - S((/7, Dv)) : ■i/’i(|qn|)D(v„ - v)d(x,o-) —)■ 0 asn^oo. 

jQt 

For t? G (0,1) we obtain 


J (^r/£»(S((/?„,Dv„) - S((^,Dv)) : D(v„ - v)^ d(x,a) 

= X{|q„|>i} Dv„) - S(<y9, Dv)) : D(v„ - v)^ d(x, a) 

+ / X{|q„|<l} (^r]Q{S{(p, Dv„) - s((^, Dv)) : D(v„ - v)^ d(x, a) 


'Qt 

= : {A) + {B). 


1 - 1 ? 


By (14.7511 and (|4.5ip there holds 

(A) < G (0,T) X spt(7/) : |qn(t,x)| > l|) 

X ^ y Dv„) - S((/?, Dv)) : D(v„ - v) d(3:, u)^ 


< c 


P dx dcr I 


1 - 1 ? 


'o 4spt(T;) 




0 as n —)■ oo, 


where we took into account Holder’s inequality. Since (B) also vanishes for n —)■ oo by 
(j4.83p . we finally have shown 

/ (,,{S(^„,Dv„)-SfeDv)):D(v„-v))" d(x, a) —)• 0 as n —)• oo, 

for all 1 ? G (0,1). We have by assumption (A2) choosing i? > ^ 

/ r/£)(S(v9„,Dv„) - S((/?, Dv„)) : D(v„ - v) d(x,(T) 

JQt 

<c [ IvJn - + |DVn|’^^^“^^)(|DVn|’^ + iDvl’^) d(x,fj) 

JQt 


< c 


' Qt 


Wn d(x,t) 


1 - 1 ? 


' Qt 


(1 + |Dv„|^ + iDvl^’) d(x, a) 


0 , 


where we have used (|4.52ll and (I4.5ip i . This finally implies 

(^7 ?£>(S((^,Dv„) - S((/?,Dv)) : D(v„ - v)^ d(x, a) —> 0 as n —oo, 

for all 1 ? G (^,l) and a.e. t G (0,T). The monotonicity of S supposed in (A2) implies 
that Dv„ —7- Dv a.e. as r] is arbitrary and g is strictly positive. This justifies the limit 
procedure in the energy integral, e.g. S = S((/9, Dv) is shown and the proof of Theorem [T] 
is therefore complete. 
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A Appendix: Pressure Decomposition 


The next theorem is in the spirit of | 33 l Theorem 2 . 6 ] but without the condition divu = 0 . 
Theorem 7 . Let u e C^([ 0 , T]; 1 < q < oo, and let H £ such that 

/ u ■ dtT] d{x,t) = / a :VT]d{x,t) 

j Qrp j Qrjp 

for all r] £ C'^(Qt)- Then there are integrable functions vr/i and ttq such that 


/ u ■ dt'nd{x,t) = / H:V77d(x,t)+ / • 

J Qj^ J Qj- J Qrj- 


ttq ■ div r] d{x,t) 


I 

J Qt 


TTh ■ dtdivr]d{x,t) 


for all T] £ C^{Qt). 

a) There holds TTh £ (^^([O, T], L™™7r/i(0) = 0 and Avr/i = — div(u — uq). 

b) We have for any 1 < r <2 


( 1 . 84 ) lk/i||L°°(0,T;L’'(n)) < c(r, g) (^||H||ir(gy) + ||u||j;^oo(o^T;L 2 (n)) j > 

( 1 - 85 ) ho\\L<J{QT) < c(g) ||H|| L<t{QT) 

for some c{r,q),c{q) > 0 independent o/u, H, tto. 

Corollary 8. Let the assumptions of Theorem be satisfied. Assume in addition that 
H = Hi + H2 where Hi £ and H2 £ for some 1 < (71,92 < 00. 

a) Then we have tto = tti + 7 r2, where 


(1-86) IkillLn(QT) < c{qi,q2)\\ili\\Lii{QT)^ 

( 1 - 87 ) \\'^2\\li2{Qj.) < 0(91,92) ||H2||l92(Qj,) 

for some constant 0(91,92) > 0. 

h) //VHj £ L'1^{Qt) for some i £ { 1 , 2 }, then we have for all LI' LI 


( 1 . 88 ) 


V7ri||ig2((0,T)xC2') < c(9i, 92) ||VHj||i;^9i(Qj,) 


for z = 1, 2 and some constant 0(91,92) > 0. 

Proof of Theorem\^ Following [ 33 ] (proof of Thm. 2 . 6 ) there is some tt £ (^^([O, Tj; L^( 0 )), 
where r = min{2,9}, with 7r(0) = 0 such that 



TT div r] d{x, t) 



(u-u(O)) 


• ?7d(a:,t) + 



Hdcj : Vr]d{x,t) 


for all 77 £ C'^(Qt)'^ and 


I'^IIl°°( 0,T;L’-(O)) < C (^||H||ir(gj,) + ||u||ioo(o,T;L2(o))^ 
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for some c > 0. Let £ be the solution operator to the bi-Laplace equation with zero 
boundary values for the solution and its gradient. Then we have that £ extends to a 
bounded linear operator 

(1.89) £ : 


for all q G (1, oo), see [30]. We decompose tt = tto -|- vr/i, where 

7ro(t) = A£F(t), {F{t),rj) = f f H(x, u) da : V^r/(x) dx for all rj ^ (Q) 

Jn Jo 

and TT/i = TT — TTo. We have 7ro(0) = 0 and hence 7r/i(0) = 0. Due to (I1.89P we have (|1.84p . 
Moreover, there holds 



for all rj G C^{Q) and so Att^ = — div(u 


u 


u(0)) • Vr]d{x,t) 


J Qrp 

-u(0)) a.e. We set TTo = dtTTQ such that 


TTo A ?7 dx = / H : dx 


for all r] G C^(D) and a.e. in t. We obtain (I1.85P as a consequence of (ll.89p . □ 

Proof of Corollary We recall from the proof of Theorem [7] that vro = ACF. So we set 


VTi = A£Fi, 7r2 = A£F2, 


where Fi,F 2 are defined analoguously to F. The claim follows from the continuity prop¬ 
erties of the operator A£ (which follow from (11.891) 1 and local regularity theory for the 
bi-Laplace equation. □ 


B Appendix: Evolution Equations for Monotone Operators 

In the following we will recall some results and definitions from the theory of monotone 
operators, subgradients and an associated evolution equation. Furthermore, we prove a 
characterization of the operator A from the proof of Lemma HI as a subgradient. For an 
introduction to the theory of monotone operators we refer to Brezis m and Showalter m- 
In the following let iL be a real-valued and separable Hilbert space. Recall that A: H ^ 
V{H) is a monotone operator if 

{w — z,x — y)H > 0 for all w G A{x)^ z G A{y). 

and V{A) = {x ^ H ■. A{x) / 0}. If (p: H ^ RU {-|-oo} is a convex function, then 
dom((/?) = {x G H : ip{x) < oo}. Moreover, ip is called proper if dom(< 7 ?) / 0. The 
subgradient dnP'- H —)■ V{H) (with respect to H) is defined by u; G dnpix) if and only if 

<f(0 > p{x) + {w,f, — x)h for all G H. 

dup is a monotone operator on H. If ip is lower semi-continuous, then dip is maximal 
monotone, cf. [131 Exemple 2.3.4]. 
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Theorem 9. Let Hq,Hi be real-valued, separable Hilbert spaces such that Hi ^ Hq 
densely. Moreover, let cp: Hq {+ 00 } be a proper, convex and lower semi-continuous 

functional such that ip = ipi -\- ip 2 , where ip 2 > is convex and lower semi-continuous, 
dom(/?i = Hi, and is a bounded, coercive, quadratic form on Hi. Set A = dny^- 

Furthermore, assume that B: [0, T] x Hi ^ Hq is measurable in t & [Oj^] Lipschitz 
continuous in v ^ Hi satisfying 


\\B{t,vi) - B{t,V 2 )\\Ho < M{t)\\vi - V 2 \\hi for a.e. t e [0,r], 

for all vi,V 2 G Hq, and for some M G L^(0,T). Then for every uq G dom((/j) and f G 
LP‘{D,T;Hq) there is a unique u G T ; Hq) n -L°°(0, T ; Hi) with u{t) G 'T>{A) for a.e. 

t > 0 solving 

nil 

(2.90) —(i) + ^(ti(i)) 3 B{t,u{t))-\-f{t) for a.a. t ^ {0,T), 

(2.91) m( 0) = Mo- 


Moreover, ip{u) G L°°(0, T). 

We refer to [H Theorem 4] for the proof. 

In order to apply the latter theorem we will use: 


Lemma 10. Let Eq be as in Section\^ Then the subgradient 81^2 Eq of Eq is single-valued 


and we have 


(0) 


V{dL 2 ^Eo) = {m G H^{n) n i/fo)(0) : dj^iplan = 0}, 
dif^^Eoiip) = -Aip + PofoiT) for all ip G V{dL 2 ^^Eo). 
Proof. The inclusion 


{m G H\n) n /7j))(0) : djpTldn = 0} C P(a^2^^i?o) 

and —Aip + f^i^ip) G 8^2 E(i{ip) can be shown in a straightforward manner using the defi¬ 
nition of the subgradient. Conversely, if m; G 81 ^ 2 ^ ^Eo{ip) for some ip G 'D{ 8 j ^2 ^Eq), then 

J w{r] — ip) dx < EQ{r]) — Eii{ip) for all rj G dom(illo) = -f^(o)(^)- 

Using this inequality with rj = ip -\-tfj ior if ^ t > 0 arbitrary, dividing by t, and 

passing to the limit t —)• 0+, one obtains 


/ 

Jo. 


w'lfdx^ / Vip ■ Vdx-\- / /o(</?)'0dx 


for all fj G ff^Q^(fl). Replacing ■0 by —0 one even obtains equality. Hence (/j is a weak 
solution of the Laplace equation with Neumann boundary conditions 


-Aip = w - Pofoiip) in 0, 

8j\fip = 0 on 8Q, 

where fQ{ip) G L‘^{Ll) due to E L®(H) and (Al). By standard elliptic regularity one obtains 
ip G H‘^{LI) and 8 j\fip\QQ = 0 in the trace sense. Hence we obtain the converse inclusion for 
'D{ 8 ji 2 ^^Eq) and 8j^2^^EQ{ip) = {—Aip -\- Po/o(‘^)}- Therefore the claim is proved. □ 
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Corollary 11. We have 

v{d^-iEo) = {u& H^{n) n ^^fip\an = o}, 

dfj-iEo{ip) = -An{-^(P + foi^)) for all G V{d^-iEQ) 

(0) -"(0) 

Proof. The statement is proved in the same way as in the proof of [5l Corollary 4.4], where 
one uses that —Ap + foip) G H^{Q) and elliptic regularity theory again. □ 
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